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Binary mixture of hard disks as a model glass former: Caging and uncaging 
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Japan 

(Dated: 30 August 2011) 

I have proposed a measure for the cage effect in glass forming systems. A binary mixture of hard disks is 
numerically studied as a model glass former. A network is constructed on the basis of the colliding pairs 
of disks. A rigidity matrix is formed from the isostatic (rigid) sub-network, corresponding to a cage. The 
determinant of the matrix changes its sign when an uncaging event occurs. Time evolution of the number of 
the uncaging events is determined numerically. I have found that there is a gap in the uncaging timescales 
between the cages involving different numbers of disks. Caging of one disk by two neighboring disks sustains 
for a longer time as compared with other cages involving more than one disk. This gap causes two-step 
relaxation of this system. 

PACS numbers: 61.43.Fs, 64.70.pv, 64.70.pm 
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I. INTRODUCTION 

When a glass-forming material is cooled, the motion of 
the constituents becomes extremely slow. Likewise, the 
mobility of a granular material decreases as the grains are 
tightly packed. Although the glass forming materials are 
thermal and the granular material is athermal, many fea- 
tures similar to that of glass-forming materials have been 
found in a granular system^. This similarity suggests a 
common mechanism in these systems. The motion of the 
constituents is restricted due to their neighbors. This is 
called 'cage effect—. Clearly, the cage effect is of geo- 
metrical origin, regardless of whether a material is ther- 
mal or athermal. If a constituent cannot escape from a 
cage after cooling or packing has taken place, the system 
freezes. This is observed as glass transition or jamming 
transition for granular materials. The cage effect plays a 
critical role in these transitions. 

The cage effect has been studied using various proce- 
dures. A plateau in the time evolution of the mean square 
displacement has been found^i between the initial ballis- 
tic regime and the later diffusive regime. A constituent 
rattles inside a cage during the plateau. The intermedi- 
ate scattering function displays a two-step relaxation in 
accordance with the mean square displacement^. Cage 
correlation function^ defined by the list of neighbor's suc- 
cessfully reproduces the temperature dependence of a dif- 
fusion constant. Three point velocity correlation function 
captures the correlation of directions of motion between 
different times, and shows a backward motion of con- 
stituents caused by a cage-. Non-Gaussian parameter 
reveals that the velocity distribution of the constituents 
deviates upwardly at a high velocity range from a sim- 
ple Gaussian distribution, corresponding to a correlated 
motion induced by a cage^. The four point correlation 
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function 7 shows the correlations between the relaxation 
of different constituents. The jump of a constituent^ can 
be divided into two types: reversible and irreversible. 
The ratio of irreversible to reversible jumps increases 
with temperature, corresponding to an increase in the 
escape from a cage. 

All of these functions or measures deal with the mo- 
tion of the caged constituents. The effects of a cage are 
discussed indirectly through the motion of the caged con- 
stituents. There is no measure directly dealing with the 
caging constituents adjoining to the caged constituents. 
I have proposed a measure for the cage effect based on 
the determinant of the rigidity matrb*2. The concept of 
rigidity has been applied to network glasses 1 ', and fluid- 
solid transition^. 

I have selected a mixture of hard disks, which has al- 
ready been investigated by many authors^—, as a model 
glass former. The disks collide with each other elastically, 
without energy dissipation. As the packing fraction cf> 
(fractional area occupied by disks) increases , a two-step 
relaxation of the density fluctuation 5 is observed in this 
system. 

In the next section, I define a local network of inter- 
acting disks, characterizing a cage. I introduce a param- 
eter which quantifies the caging in Section 3. Numerical 
methods are explained in Section 4. Numerical results 
are presented in Section 5. Section 6 gives discussions 
and conclusions. 



II. LOCALLY RIGID FRAMEWORK 

Figure 1 shows two examples of the cage effect. The red 
disk caged between two black disks [Fig. 1(a)] cannot pass 
through between the two disks to the other side unless the 
gap between the two disks increases up to the diameter 
of the red disk. Likewise, a triangle of three disks con- 
strained by the three surrounding disks [Fig. 1 (b)] cannot 
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FIG. 1. (color online). Schematics of the cage effect, (a) 
Red disk is caged between two black disks. The red disk 
cannot pass through between the two disks (from the left to 
the right) unless the gap between the black disks is larger than 
the diameter of the red disk (middle). Fictitious networks 
representing contacts between disks are shown by black lines, 
(b) A cluster of three blue disks is caged by three black disks. 



rotate unless the top disk shifts upwards such that the 
two triangles, one of which is formed by the centers of 
the three caged disks, and the other consists of the cen- 
ters of the surrounding disks, are perspective from a line 
(the extensions of three pairs of corresponding sides of 
the two triangles meet in collinear points^). 

If we apply an external force to the red disk to make it 
pass through between the two caging disks or to the blue 
disks to rotate the triangle, force carrying networks are 
formed as shown by the black bonds in Fig. 1 provided 
that the positions of the end disks are fixed. If we place 
hinges at the centers of the disks connected by rigid bonds 
rotating freely around the hinges, the networks in Fig. 1 
are rigid and isostatioi^, such that any removal of a bond 
results in a network that is not rigid. If the number of the 
caged disk is Adisk , the number of columns of the matrix 



is 2iVdi s k (there are two bonds in Fig. 1(a) (Adisk = 1) and 
six bonds in Fig. 1(b) (A^isk = 3)). This is a property of 
an isostatic network. 

In the disk system, the combination of the collid- 
ing pair of the caged disks does not change for a short 
timescale. From the list of colliding pairs, we can con- 
struct a network where each bond represents the collid- 
ing disk pair. If a disk (or a set of disks) is caged, we 
can find a set of bonds corresponding to those shown in 
Fig. 1. This is because the condition of uncaging shown 
in Fig. 1 is the same even if the disks move and collide 
with each other. Therefore, we can define an isostatic 
sub-network in the network clS Si C9. ge. This local sub- 
network of bonds is called the 'locally rigid framework 
(LRF)'. An LRF inhibits the motion of a disk (or a set 
of disks) in one particular direction. If all LRFs pertain- 
ing to a disk are effective, the disk can only rattle inside 
the cage. It should be noted that this network is not 
a real contact network, but a fictitious one because the 
disks do not continue to contact but they collide. More- 
over, an LRF may be broken as the system evolves and 
distances between the disks increase. 



III. CAGING PARAMETER 

An uncaging event can be detected by the rigidity 
matri^. We can construct the rigidity matrix from the 
LRF such that each row represents a bond and each col- 
umn corresponds to the coordinates of caged disks. An 
LRF consists of bonds inside caged disks in addition to 
bonds connecting between a caged disk and a neighboring 
caging disk. There are two bonds connecting a caged disk 
(red) and caging disks (black) in Fig. 1(a). For Fig. 1(b), 
there are three bonds inside the caged disks (blue) and 
three bonds connecting the caged disks and caging disks 
(black). 

If there are bonds between the caged disks 1-2 and 2 
- 3, the rigidity matrix is written as^ 
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where Xi and yi are the positions of the i-th disk, and 
dij = y/ {xi — Xj) 2 + (yi — yj) 2 is the distance between 
the centers of disks i and j. It should be noted that the 
matrix is not symmetrical with exchanging i and j. If i- 
th disk is caged by j-th caging disk, the matrix contains 
(xi — Xj)/dij but does not contain (xj — Xi)/dij. 

The rigidity matrix for an isostatic network is a square 
matrix (the number of degrees of freedom of the caged 



disks 2A r disk and the number of bonds are the same in an 
isostatic network 9 ) and has a determinant except for a 
special position of the disks 9 - when uncaging occurs. 

When the red disk passes between the other two disks 
in Fig. 1(a), the determinant of the LRF becomes zero 
because the two bonds are in parallel when the three cen- 
ters of the disks are collinear. The determinant changes 
its sign at this point. Likewise, the sign of the determi- 
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nant in the right panel is different from that in the left 
panel in Fig. 1(b). 

Using the determinant of the LRF with A^isk caged 
disks, we can define the caging parameter XN disk , which 
is a measure of the effectiveness of cages. It is defined as 



where, £jv diak i s the number of LRFs investigated which 
have -ZVdisk caged disks, and nAr diBk (t) is the number 
of LRFs whose determinants change their sign among 
SAr disk . It should be noted that the number IIjv dIsk (i) 
is increased by one when a determinant changes its sign, 
only if an LRF has not been broken such that the lengths 
of all the bonds d,j are less than f(Ri + Rj), where Ri is 
the radius of the ith disk. The factor / = 1.5 is adopted 
here. 



IV. NUMERICAL METHODS 




FIG. 2. (color online). An example of LRFs. Disks and bonds 
are chosen in a circular region in the computational area. Sets 
of bonds denoted by different colors are identified as LRFs on 
the basis of the rigidity of the local sets of bonds. An LRF of 
6 bonds is denoted by the thick pink lines. The thick black 
line is the residual bond after decomposition 



To determine the evolution of Xjv diak (i)j the positions 
of disks should be known. The evolution of the posi- 
tions of hard disks can be efficiently computed by an 
event-driven simulation 14 , in which the collision times 
between disks and velocities after collisions are sequen- 
tially solved. Disks having two different radii, 1.0 and 
1.4 in dimensionless units, are mixed in equal propor- 
tions. This size ratio has been chosen to safely avoid 
crystallization 1 ^. The total number of disks is 1024. A 
periodic boundary condition is applied. The packing 
fraction of the disks <f> is 0.8. The initial velocity dis- 
tribution is Gaussian with unit variance in dimensionless 
units. I have defined the colliding pair of the disks in 
a particular time window of 200 in the normalized unit. 
The fictitious network is constructed from the list of col- 
liding disks in the time window. 

In order to construct an LRF, a circular region is first 
placed randomly on the fictitious network. The disks 
and the connecting bonds inside the circle are chosen. 
By varying the radius of the circle, the number of disks 
inside the circle is changed between three and hundred, 
which is much smaller than the system size (1024 disks). 
The number of chosen bonds is usually more than that 
required for an isostatic network. This is because the 
disks move and have a chance to collide with more disks 
than when they are tightly packed. 

Overconstraining bonds (their removal does not affect 
the rigidity of the network) are randomly removed until 
the isostatic condition 1 ^ is fulfilled. An example of an 
isostatic network is shown in Fig. 2. If a bond is further 
removed from the isostatic network, some fraction of the 
bonds becomes mobile. Correspondingly, a disk (or a 
set of disks) also becomes mobile (rotates around a disk 
center). For example, if we remove a blue bond at the 
bottom of Fig. 2, another blue bond becomes mobile and 
the disk connected by this blue bond can rotate around. 
We can find a minimum set of mobile bonds among 2 A^isk 



possible bond removals. The mobile bonds can be found 
by the pebble game algorithm 1 ^. 

In addition to the removed bond, the induced mobile 
bonds are further removed from the isostatic network. 
In the example above, we remove two blue bonds from 
the network. An LRF is composed of this local mobile 
network plus the removed bond. The remaining net- 
work, after the removal, is still isostatic. During each 
removal of bonds, the mobile disk (caged disk) is also re- 
moved accordingly. The disks adjacent to the mobile disk 
connected through the removed bonds correspond to the 
caging disks. This removal process is repeated. Finally, 
the isostatic network is uniquely decomposed into the sets 
of LRFs plus one bond (thick black bond in Fig. 2). In 
Fig. 2, there are 17 LRFs with 2 bonds (Adisk = 1) and 1 
LRF with 6 bonds (Adisk = 3) shown by pink lines. The 
pattern of LRFs shown in Fig. 2 changes if the position 
of the randomly chosen circular region is shifted or the 
random sequence of the removal of the overconstraining 
bonds is changed. 

According to this procedure, the LRF of A^isk = 2 is 
impossible because we can always find an immobile bond 
for a particular removal of a bond. An LRF of Adisk = 2 
splits into two LRFs of A^isk = 1. 

V. NUMERICAL RESULTS 

A. Two-step relaxation 

Black curves in Fig. 3 display the evolution of XN diBiL - 
From the right to the left of the black curves, the num- 
ber of caged disk Adisk is 1, 3, 5, and 10. As Adisk in- 
creases, XAr diBk decreases faster. All %s start to decrease 
at t = 10~ 2 , when the intermediate scattering function 
F qmax (t) starts the first relaxation. F qniax (t) is given by 
(Si exp(g max (a; i (t) — sGj(0)))), where Xi is the position 



Time N 
o I , , , , llisJ 

10 2 10° 10 2 10 4 10 6 

Time 

FIG. 3. (color online). Decrease in caging parameter XJv diBk (t) 
for various iVdisk at <f> — 0.8. From right to left, the black 
(solid) curves correspond to iVdisk = 1, 3, 5, and 10. The 
red (upper dotted) curve represents the subset of iVdisk — 1 
in which a bond exists between the two caging disks when 
the bonds are placed according to the colliding pairs. The 
green (lower dotted) curve corresponds to the other group 
in which no bond is placed between the caging disks. The 
orange (dotted-dashed) curve is the same as the red (upper 
dotted) curve but njv diak (t) in Eq. (2) is multiplied by a factor 
of 20. The blue (dashed) curve is the intermediate scattering 
function (self-part) F qjnsjc (t) at the wave number q ma x where 
structure factor is maximum, [inset: mean square displace- 
ment] 



of the i-th disk, the angular bracket is the ensemble av- 
erage, and g max is the wave number at which the struc- 
ture factor is the maximum. This function measures the 
relaxation of the structure with the nearest-neighboring 
length scale. This first period corresponds to the ballistic 
motion of the disks before a collision as can be seen in 
the inset of Fig. 3 where the mean square displacement, 
defined by {[x{(t) - cc^O)] 2 }, increases as t 2 . The LRFs 
uncaged during this stage are essentially ineffective, or 
very weak. 

F, mix (t) reaches a plateau when t ~ 4. All %s continue 
to decrease after F 9majt (t) reaches the plateau. We cannot 
find any coincidence between \ s (solid curves) and the 
second relaxation of F qin ^(t). 

The decrease in X3, xs> and Xio ceases when t ~ 10 2 . 
All LRFs belonging to these %s are uncaged or broken by 
this time. 

The decrease in xi (the rightmost black curve) is 
clearly distinct from the others. There is a gap in the 
uncaging times between xi and X3- The uncaging time 
(the time when XN disk (t) = 0.9) for xi is larger than 
that for X3 by a factor of 1.8 x 10 3 (Fig.4(a)). It should 
be noted that there is no particular behavior in xi at 
t ~ 10 2 , when i 7 'g max (t) starts its second relaxation. 

This gap is emphasized if we divide the LRFs of 
Adisk = 1 into two groups. In one group Xi,lo, there 
is a bond between the two caging disks, when the bond 



•1 1 , , , 

5 10 15 20 

AAdisk 




FIG. 4. (color online). Gaps between cages, (a) Uncaging 
time as a function of iVdisk. (b) The distribution of the ab- 
solute value of the eigenvalue e of the rigidity matrix. The 
red curve stands for iVdisk = 1 with a bond between the two 
caging disks, and the green curve for iVdisk = 1 without the 
bond. The two black curves are for iVdisk = 5 and iVdisk = 10. 



is placed on the basis of the colliding pairs. In the other 
group Xi.wk, no bond is placed between the caging disks 
(schematics are shown in Fig. 4(b)). The decrease in 
Xi,lo a nd Xi,wk are shown as the red (upper dotted) 
and green (lower dotted) curves respectively, in Fig. 3. 
The green curve starts to decrease at the same period 
to other xs and continues to decrease through £ ~ 10 2 
without any notable change. 

On the other hand, Xi.lo starts to decrease at t ~ 10 2 . 
Interestingly, this period coincides with the onset of the 
second relaxation of F qmax (t). There is no special sig- 
nature at t = 10 2 in other xs. The amount of the de- 
crease in xi.lo is small, because more than 90% of the 
LRFs investigated are broken before uncaging. The or- 
ange (dotted-dashed) curve is the same as the red curve, 
but IlAr disk (i) in Eq. (2) is multiplied by a factor of 20 
to make a comparison with F Qmgx (t). A good coincidence 
can be seen between the second relaxation of F qmax (t) 
and xi,lo- 
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B. Gap between cages 

Figure 4(a) shows the uncaging time (the time when 
XN disk (t) = 0.9) as a function of Ad; s k- It can be seen that 
the uncaging time of LRFs with Adisk = 1 is significantly 
larger than those of other x s - The uncaging time of other 
%s decreases with Ad; s k- A determinant M of an LRF can 
be factorized such that M = HS, where H is composed 
of the bonds of an isostatic subgraph (a triangle) in an 
LRF, and S depends on the pattern of the rest of the 
LRF—. A determinant of an LRF can be decomposed 
into factors corresponding to triangles and to connecting 
bonds between triangles. A factor for a triangle cannot 
change its sign because the LRF is broken before the 
change. Therefore, the change in sign of a determinant 
for iVdisk > 3 comes from the factors of connecting bonds 
between triangles. Because the number of such factors 
contained in M increases as Adisk, the uncaging time 
decreases as Ad; s k- 

The distribution /(e) of the absolute value e of a com- 
plex eigenvalue of the rigidity matrix when i = is shown 
in Fig. 4(b). /(e) is defined such that f(e)de is the frac- 
tion of the eigenvalues between e and e + de and that 
/ fde = 1. The eigenvalues evolve with time followed by 
the evolution of the disk position. When one of the eigen- 
values of a rigidity matrix goes to zero, the determinant 
changes its sign. Thus, the eigenvalue e measures the 
difficulty of the uncaging. As e increases, the displace- 
ments of the disks required for uncaging increases. An 
eigenvalue corresponds to an eigenfrequency if there is a 
harmonic potential between the disks. It can be seen that 
there is a substantial fraction of low frequency modes for 
Adisk > 3 in addition to Adisk = 1 without a bond. On 
the other hand, the distribution /(e) for Adisk = 1 with 
a bond falls to zero at e = 0.3. This difference is the 
origin of the gap in uncaging times. An LRF with small 
e uncages faster. However, an LRF consisting of tightly 
packed disks sustains for a long time. 

VI. DISCUSSION AND CONCLUSION 

These observations suggest that the uncaging of LRFs 
of Adisk = 1 with a bond (xi.lo) requires the uncaging 
of other LRFs. Uncaging of LRFs with Ad; s k > 3 and 
Adisk = 1 without a bond proceeds during the plateau of 
the intermediate scattering function. The LRFs relaxed 
in the first stage flip over again between the neighboring 
configurations as shown in the left and right panels in 
Fig. 1. The system gradually changes its configuration 
by the flipping of the weak LRFs relaxed in the first stage. 
The plateau shows that the disk system requires time to 
find a configuration to open a path between two disks 
packed initially with an order. 

Displacements of disks due to the uncaging of these 
LRFs are small, and can be seen as a plateau of the 
mean square displacement when t ~ 10, as shown in 
the inset of Fig. 3. This period corresponds to /3 relax- 



ation of glass forming materials^. Disk positions gradu- 
ally change in this period and the distance between two 
caging disks in LRFs with Adisk = 1 with a bond, en- 
larges. The uncaging of LRFs with Adisk = 1 with a 
bond triggers a large structural change and promotes the 
second relaxation of F qmax (t). This period corresponds 
to a relaxation of glass forming materials^. The require- 
ment of uncaging of other types of LRFs for LRFs with 
Adisk = 1 with a bond is the situation of "constrained 
dynamics"— which explains the slow relaxation of glass- 
forming systems. According to this study, uncaging of 
LRFs Adisk = 1 with a bond cannot take place by it- 
self. Through the uncaging of LRFs with Adisk > 3 and 
Adisk = 1 without a bond, the arrangement of the disks 
changes so that uncaging of LRFs Ad; s k = 1 with a bond 
takes place. 

An LRF with Ad; s k = 1 with a bond is the collection of 
locally ordered packing of three disks because the three 
disks collide with each other when the network is deter- 
mined. This local order cannot propagate throughout 
because of the poly-dispersity of the disks. This situa- 
tion is called the geometrical frustration-^. Because of 
the lack of the long-range order, the geometrical frus- 
tration inevitably leads to the formation of weak LRFs, 
which are uncaged first. 

In conclusion, I have proposed a measure for the cage 
effect of a binary mixture of disks. I have found the 
gap between cages with different numbers of caged disks. 
Cages with one disk relax slower than those involving 
more than one disk. The former cages correspond to 
the locally ordered packing of disks. The latter cages 
represent disordered packing of disks. This coexistence 
of locally ordered and disordered disks leads to the two 
step relaxation. 
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